NULL CONTROLLABILITY OF ONE DIMENSIONAL 
PARABOLIC EQUATIONS 

G. ALESSANDRINI AND L. ESCAURIAZA 

Abstract. We prove the interior null-controllability of one— dimensional par- 
abolic equations with time independent measurable coefficients. 



1. Introduction 

Let us consider the following one-dimensional heat equation with variable coef- 
ficients 
(1.1) 

d x (a(x)d x z) + b(x)d x z + c(x)z — p(x)d t z — fxu, 0<x<l, < t < T , 
< z(0,t) = z(l,t) = 0, < t < T , 

z(x,0) =z , < x < 1 . 

Here, z(x,t) is the state and f{x,t) is an interior control that acts on the system 
over the open set u) C (0,1). The coefficients a, b, c and p are assumed to be 
measurable, bounded and for some K > 1, 

(1.2) R- 1 < p{x) < K, R- 1 < a(x) < K, \b(x)\ + \c{x)\ < K , a.e. in [0, 1] . 

For any given zq in L 2 (Q, 1) and / in L 2 (uo x [0,T]), there is only one solution 
z to (O) in C([0,T];L 2 (0,1)) n L 2 (0, T; Hq(0, 1)). The goal of this paper is to 
analyze the interior null-controllability of (|l.l(l . Specifically, we want to solve the 
following problem: 

Given T > and zq in L 2 (0, 1), to find f in L 2 (uj x [0, T}) such that the corre- 
sponding solution z to (|l.lfl satisfies, z( ■ , T) = in (0, 1). 

In h is shown that the system (|1.1|) is null-controllable at any positive time, 
when the coefficients a and p are Lipschitz in [0, 1]. In this reference, the proof of 
null-controllability is based on an appropriate observability inequality for the ad- 
joint system and it is implied by a global Carleman estimate. When the coefficients 
are smooth, the observability inequality can be proved introducing Fourier series 
and using high frequency asymptotic formulae for the eigenvalues of the correspond- 
ing Sturm-Liouville problem and classical results on the sums of real exponentials, 
see jEj. In [0] adopting the approach introduced by D.L. Russel in ^5] (the null 
controllability of the wave operator at large times implies the null controllability of 
the heat equation at all times) it is shown that the system 1)1.111 is null controllable, 
when a and p have bounded variation in [0, 1]. 

The main result in this paper is the following. 
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Theorem 1. Assume that the coefficients a, b, c and p are bounded measurable 
and satisfy ()1.2|) . Then, Ijl.lJ) is null- controllable at time T , for all T > and with 
controls f in L 2 (0, T; H^(uj)). 

To prove this result we proceed in the the following way. First, a change of 
variables shows that the internal controllability of the system is equivalent to 
the same question for a system 

!d 2 x z - p{x)d t z = fxu, < x < 1 , 0<t<T , 
z(0,i) =z{l,t) = 0, 0<t<T, 
z(x,0) =z Q , < x < 1 , 

where p is a new measurable function satisfying l|1.2|l . for some new constant only 
depending on K, which we shall continue to denote by K. 

Then, if we denote by S the inradius of the open set lj, that is 

(1.4) 5 = sup{r > 0| 3£ G lj , g-r,(|r)cw), 

and if e\, • ■ ■ , e n . . . and < \ 2 < \ 2 < ■ ■ ■ < \ 2 n . . . are respectively the 
eigenfunctions and eigenvalues of the Stiirm-Liouville problem 



(1.5) 



we prove the following: 



e" + p(x)X 2 e = , < x < 1 
e(0) = e(l) = 0, 



Theorem 2. Assume that the coefficient p satisfies (11.211 . Then, there is a constant 
N , which depends on K and on 5 such that the inequality 

(1.6) Y a 2 < Ne N » [ \ V a k e k \ 2 dx , 



A fc <M JX \h<H 

holds for all fi > 1 anc? aZ/ sequences {a k }- 

The analog of this inequality for the eigenfunctions of the Laplace-Bcltrami op- 
erator on a compact and smooth Riemannian manifold with a possibly nonempty 
boundary was proved in |14j . There, G. Lebeau and L. Robbiano showed that it 
implies the interior null-controllability of the heat equation over the manifold by 
giving an explicit construction of the control function / (See ^3 §5] for a more 
simplified presentation). 

The arguments in |14) show that the same iterative method of construction of 
the control function / given in 14 works for the system 1|1.3|) . when Theorem [21 
holds. Thus, Theorem ^ follows from Theorem 

To prove Theorem El we start by following the arguments in |14j . In particular, 
given p > 1 and a sequence of real numbers oi, a^, ■ ■ ■ , a n . . . , we set 

u(x,y) = y a k e k (x)cosh(X k y) . 

This function satisfies 

' d\u + d y (p(x)dyu) = 0, D<Kl, yeR, 



(1.7) 



u(0,y) = u(l,y) = 0, yet, 
d y u(x,0) = 0, < x < 1 



NULL-CONTROLLABILITY OF 1-D PARABOLIC EQUATIONS 



3 



and the proof of Theorem is a consequence of a quantification of the following 
qualitative result of unique continuation from the boundary: 

Assume that u satisfies i|1.7fl and u{x,§) = ; when x is in u> C (0, l).Then, 
u = in [0, 1] x R. 

In |14j . the one dimensional interval [0, 1] is replaced by a compact and smooth 
manifold M, <9 2 by the corresponding Laplace-Beltrami operator on M and the 
authors work out the quantification of a similar qualitative property of boundary 
unique continuation for the elliptic operator, A + dy, where A is the Laplace- 
Beltrami operator on M. To carry out this quantification they use two Carleman 
inequalities. Those methods require that the elliptic operator involved has Lipschitz 
second order coefficients and so, they can not be applied to the elliptic operator in 
(II. 7|) . which has measurable coefficients. 

On the other hand, if a is a 2 x 2 symmetric and measurable matrix in the plane 
verifying the ellipticity condition 

(1.8) K- x \^ <a(x,y)Z-£<K\Z\ 2 , when (x,y) and £ 6 K 2 , 
the weak solutions of the equation 

(1.9) V ■(a(x,y)Vu) = , 

satisfy the strong unique continuation property: 

If a W{£ -solution of (J3~HJ| on a connected open set f2 has a zero of infinite order 
at an interior point, then it must be zero. 

See |21 . This qualitative result of strong unique continuation for uniformly elliptic 
equations in two independent variables is based on the connection between the 
solutions of these equations and the theory of quasiregular mappings and on 
the so-called Ahlfors-Bers representation |TJ of such mappings. Here, we describe 
some quantifications of this qualitative result and apply them to prove the null- 
controllability property. In particular, a "Hadamard's three circle theorem" , 
Proposition ^ an d a "doubling" type property, Proposition [21 adapted to the 
solutions of H1.9(l . 

In section [21 we recall the results we need from the theory of quasiregular map- 
pings and prove the adapted Hadamard's three circle theorem and doubling prop- 
erty. In section |21 we show how to apply them to prove Theorem [21 also using an 
estimate of continuation from Cauchy data for solutions of (|1.9fl Lemma ^ which 
we adapt from [3]. It may be worth noting that the approach used for the proof 
of Lemma Q is based on a variation on the classical principle of majorization by 
harmonic measure, |2()l Chapter VIII, §1, p. 301], which in turn has its roots in 
arguments due to Carleman p. 3-4]. 

2. Quantitative estimates of unique continuation with discontinuous 

coefficients 

Throughout the paper, z — x + iy, f2 is a simply connected open set in the plane, 
B r a circle of radius r centered at the origin, and 

= \ (0 X f + idyf) ,d g f=\ (d X f - idyf) . 

We shall denote by C constants only depending on K, whereas by N we shall denote 
constants only depending on K and 5. 



4 



G. ALESSANDRINI AND L. ESCAURIAZA 



When u 6 W^„ c (fi) is a weak solution to H1.9(l . and cr satisfies ()1.8|l we can 
associate in a natural fashion, which generalizes the harmonic conjugate, a new 
function, the so called stream function v, which satisfies 

(2.1) Vv = JaVu 
almost everywhere in f2 and is a weak solution to 

(2.2) v(g^W)=0, infi. 

Here J denotes the matrix representing a 90° rotation in the plane 

-1 



J 



1 



Moreover, letting / = u + iv, we have / £ W. ' and satisfies 

(2.3) dzf = fid z f + vd z j ', almost everywhere in ft , 

where the complex valued functions fj, and v can be explicitly expressed in terms 
of a , see 0, and verify 

(2.4) + \u\ < jfcj < 1, almost everywhere in f2 

That is, / is a i^-quasiregular mapping. 

To give an idea of why these results hold, observe that the vector field 

is, in the weak sense, curl- free in Q. To verify that v is a Wi ' -solution of l|2.2|l . 
observe that, from l|2.1|l . one obtains that the vector field 

-£-Vv = J*- 1 J t Vv 

is, in the weak sense, divergence-free in Q. 

By the Ahlfors-Bers representation (see also [5] and 0] Chapter II. 6, pp. 
258-259]), any K-quasiregular mapping / in B\ can be written as 

f = Fo X , 

where F is holomorphic in Bi and £ = x( z ) 1S a -K'-quasiconformal homeomorphism 
from B\ onto B\, which verifies, x(0) = 0, x(l) = 1) 

(2.5) C- 1 \z l -z 2 \^ < \x(zi) - X(Z2)\ <C\ Zl -z 2 \ a , whenz 1 ,z 2 eB 1 

for some < a < 1 and C > 1 depending only on X. 

We now recall the Hadamard's three-circle theorem |18) . 

Theorem 3. Let F be a holomorphic function of a complex variable in the ball B r2 
and M(r) = maxg r \F\. Then, the following is valid for < n < r < r 2 , 

log ^ log 
log M(r) < -^logM(n) + — £ log M(r 2 ) 
log log : 



The meaning of this inequality is that log M{r) is a convex function of the 
variable logr. 

Let u <E W (Br) be a weak solution to l|1.9f) and let / : Br — > C be the 
associated if-quasiregular mapping. Rescaling (|2.5(l we have that f — F ox, where 
F is holomorphic in Br and x '■ Br — ► -B^ is a i-T-quasiconformal homeomorphism, 
which verifies 

(2.6) RC-%\* < \ X (z)\ <RC\^\ a 



NULL-CONTROLLABILITY OF 1-D PARABOLIC EQUATIONS 



5 



where C is the same as in ()2.5|) . 
Define 

(2.7) B r = {z G B R : \ X (z)\ < r} 
and 

(2.8) m(r) = rnax|/(z)| , when r < R . 

Then, through the change of coordinates, C, = x{ z )i the Hadamard's three circle 
theorem takes the form: the function log m(r) is a convex function of logr, 



log ^ log f- 

(2.9) logm(r) < - — ^-logm(ri) + j- — £Mog m(r 2 ) , when < r\ < r < r 2 < R 
og ri og ri 

and the sets B r , almost look like balls. In particular, 

(2.10) Br = Br , B i C B r C B R(e ,y , when r <R 

and C is the same constant appearing in l|2.6|) . Note incidentally that (|2.9|l implies 
a weak unique continuation property, that is, if m(ri) = for some small r\, then 
m(r) = for all r < R. 

On the other hand, the difference quotients of convex functions are nondecreasing 
functions of their arguments. This implies that, if / is not identically zero, 

logm(f) -logTO(|) < logm(f ) -logm(-f) < 

log § -log 5 " log f- log f 

and thus, 



77l(0 TO 



(2.11) — jjrr < —j^- , when r < R . 



m(5) 



We may prescribe that the the stream function v of u satisfies v(Q) = 0. We have 
that F = u + iv is holomorphic in the ( = ^ + ii] coordinates in Br, hence, solving 
the Cauchy-Riemann equations, 

w(C>»?)= / "c(£;,s)gLs- / u n (t,0)dt, in B R . 

JQ JO 

This formula and interior estimates for harmonic functions |llj show that in the 
^-coordinates we have, 

R 

||u|U°°(B r ) < max |F(C)| < C\\u\\ L ^ iB2r) , when r < — 

B r Z 

where C > is an absolute constant. In the z-coordinates, the last inequality reads 
as 

R 

(2.12) |M|z«>(s r ) < max\f(z)\ < C\\u\\ L ^ {B2r} , when r < — , 

and from l|2.12|l . (|2.9|l we obtain: 

Proposition 1. Let u G W,' (Br) be a weak solution to l|1.9|l and let B r , < 
r < R be the open sets introduced in (|2.7() . i/ien we have 
(2-13) 

IMU°°(B 5 ) < C ll M lli~(B r jll u lli~% r3 ) . when ri <r <r 2 < R , 6 



log^ 



'J/ ■' "~ ^r ir ■ - log — 
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And also, from [j2.11jl : 

Proposition 2. Let u and B r be as above. If u is not identically zero, then we 
have 

(2.14) ph^iEA < c NI*-(b«) , when r < fl . 

2 4 

These are respectively a Hadamard's three circle theorem and a doubling prop- 
erty adapted to the solution u through the family of "balls" B r . 

3. Proof of Theorem [I] 

First let us note that, possibly replacing z in (|l.l|l with e~ 2Kt z, we may assume 
that c is nonpositive. Introducing B(x) — f£ ds , we observe that (|l.lfl can be 
rewritten as 

e~ B Wd x (a{x)e B ^d x z) + c(x)z - p(x)d t z = f X » , in (0, 1) x (0, T] . 
The solution w to 

U~ B{X) 1 («We B(l) t) + c(x)w = , 
\w(0) = w(l) = 1 , 

verifies, < w(x) < 1 in [0, 1] and replacing z with the new dependent variable 
z = z/w, which we denote again z, we have 

e~ B ^d x (a{x)w 2 (x)e B ^d x z^ - p(x)w 2 {x)d t z = w{x)f X u, , in (0, 1) x (0, T] . 
Setting 

L= [ a- 1 (s)w- 2 (s)e- B{s) ds , y=± [ a - 1 (s)w- 2 (s) e - B ^ ds 
Jo Jo 

and writing z(x, t) — z(y,t), the new function z, which again we rename z, is a 
solution of a system of the form l|1.3fl • 

Considering the associated Sturm-Liouville problem (|1.5fl . we extend the eigen- 
functions ej, j > 1, to [—1,1] by an odd reflection in 0, similarly we extend p 
by an even reflection in 0. Next, we continue these new functions to all of R as 
periodic functions of period 2. The extended p verifies l|1.2l> . &j £ C 1,1 (M) and 
e'j + p(x)X 2 ej = 0, almost everywhere in R. 

Being the change of variable y = y(x) bi-Lipschitz, with Lipschitz constants 
which only depend on K, the open set w is transformed into a new open subset of 
(0, 1) whose inradius is comparable to 6. We continue to denote the transformed set 
and its inradius by u> and (5, respectively. Also we can assume, up to a translation 
along the real line, (—6, J)cwC (—1, 1). 

Given p, > 1 and a sequence of real numbers a\, 02, ■ ■ ■ a n , . . . , the function 



verifies 
(3.1) 



u(x,y) = ^2 a kek{x) cosh(Afcy) 



d 2 u + d y (p(x)d y u) = , in R 2 



d y u(x, 0) = , in 
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and its stream function v can be chosen so that, v £ W io ' c (K 2 ), 

\d x v = -p(x)d y u , ^ (d x {^d x vj + d%v = , in K 2 , 
[a y w = d x u , \v{x, 0) = 0, in R . 

Let / = u + iv, consider the family of "balls" B r associated to / in section [2 at 
scale R and choose R = 2 (4C) ~ , where a and C are the constants in (|2.(i[l and 
(|2.1U|I . With this choice, Br D B%. The interior bounds for subsolutions of elliptic 
equations give 

HuIU-CBjO < ilklU 2 (s 2H ) • 
These and the orthogonality of the eigenfunctions ej , j >1 imply that 

(3.2) ll»IU-( B ,) < e c, 

l|w|U«>(Sa) 

4 

An iteration of (|2.14|l and l|3.2|l give 

||«IU-(B S ) <e Ck »\\u\\ L ~ {B r) , whenfc> 1 

and from 

JOB, there is fc = k(S,K) such that, Sjl C Bs. Thus, 

2 fc 2 

(3-3) \\u\\ La , (Bl) <e N ^\\u\\ L ^ (Bs) . 

2 

The following inequality, which is an estimate on the continuation from Cauchy 
data, holds. 

Lemma 1. There are constants < 8 < 1 and C > 0, only depending on K , such 
that the inequality 

ll«IUoo C B 5 )<Cr-*|K.,0)||i a( _ riP j||«[[i- s (B4r j 

holds, when r < 1. 

The Lemma and (|3.3|) give 

HU-CBO < Ne N »\\u(- ,0)\\ L 2 { _ StS) , 

which proves Theorem 

Proof of Lemma^\ The proof is essentially contained in Theorem 4.5]. For the 
sake of completeness we summarize here the argument. 

Recalling that / = u + iv is analytic in the C-variable, we have that log |/(C)| is 
subharmonic in Br and consequently one can verify that log is a subsolution 

for an elliptic operator in divergence form E with a matrix of coefficients verifying 
l]1.8[l . For r > 0, let w be the solution to 

!Ew = 0, inB+, 
w = 1 , in (— r, r), 
u> = , in dB+ \ (— r, r) . 

On dB+ we have 

log \f\ < w\og |K- , 0)|| £ - ( _ r)r) + (1 - Ifl) log |]/|U»(B r ) 

and the maximum principle implies that the same inequality also holds in Bf . The 
Holder continuity at the boundary of w and the Harnack's inequality 11 show that 
there is 77 £ (0, 1), which only depends on K, such that w(z) > 77 in Bt . 
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Using v(0) = 0, p.ljl and interior bounds for elliptic equations [H], we have 

IMU- (Bp ) < C\\Vv\\ LHB2r) < C\\X7u\\ L2(B2r) < C\\u\\ L - (B4r) . 
These imply 

(3-5) HU°°(s 5 ) < C\\u(- ,0)\\l x{ ^ rr) \\u\\]^ iBir) . 

For every a G (0, 1], we have the interpolation inequality 

(3.6) |MU°o(-r,r) < ^(iblli^-r.rjIvlc^C-r.r) + IMU 2 (-r-.r)) , 

where (3 — 1 ^ a , C > only depends on a and |y>|c«(-r,r) denotes the standard 
C a seminorm. Next, we use the interior Holder bound for u, 

(3-7) \u\c*(b t ) <Cr- a \\u\\ La o (Bir) , 

with C > and a G (0, 1] only depending on K. Combining (|3.5|l with i|3.6|) and 
H3.7[) . we obtain the thesis with = 

The interpolation inequality l|3.6[l can be proved essentially along the same lines 
as the interpolation inequalities in [111 §6.8]. 

□ 
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